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Abstract 

In the heavy quark hmit, properties of heavy mesons and doubly heavy baryons are related by 
heavy quark-diquark symmetry. This problem is reanalyzed in the framework of Non-Relativistic 
QCD (NRQCD). We introduce a novel method for deriving Potential NRQCD (pNRQCD) La- 
grangians for composite fields from vNRQCD, which contains quarks and antiquarks as explicit 
degrees of freedom and maintains manifest power counting in the velocity via a label formalism. 
A Hubbard-Stratonovich transformation is used to eliminate four quark interactions in vNRQCD 
and then quarks and antiquarks are integrated out to get effective Lagrangians for composite fields. 
This method is used to rederive Lagrangians for the QQ and QQ sectors of pNRQCD and give a 
correct derivation of the 0(l/mQ) prediction for the hyperfine splitting of doubly heavy baryons. 
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Mesons containing a single heavy quark (Qq) and baryons containing two heavy anti- 
quarks (QQq) are related by quark-diquark symmetry. The two heavy antiquarks experience 
an attractive Coulomb force which can be computed in perturbation theory and the ground 
state of this two body system is a diquark in the 3 of color with spin 1. The mass of the 
diquark is very large and the spatial extent of the diquark is small, r ~ l/{mQv) <^ 1/Aqcd, 
thus the light degrees of freedom in the baryon orbit a static source of color in the 3 repre- 
sentation. The light degrees of freedom in the heavy meson also orbit a static source of color 
in the 3 representation, therefore the configuration of the light degrees of freedom in the two 
hadrons is identical. The diquark and quark, however, have different spin-dependent inter- 
actions which are suppressed in the heavy quark limit. Furthermore there are corrections 
coming from the composite nature of the diquark. 

In Ref. these ideas were implemented in the framework of heavy quark effective theory 
(HQET) by writing an effective Lagrangian for heavy quarks and diquarks which exhibits 
a f/(5) superflavor symmetry that permutes the five spin states of the heavy quark and 
diquark: 

C = hUDoh + iDo T + -^h^ (T-Bh- -^T^ ■ B x T + ... . (1) 

2mQ 2mQ 

Here h is the spin-| heavy quark field, T the spin-1 diquark field, and B the chromomagnetic 
field. The first two terms are the leading order Lagrangian which is manifestly U (5) symmet- 
ric. The last two terms are the heavy quark spin-symmetry breaking ©(l/mg) corrections 
which are responsible for the hyperfine splittings of the ground state heavy mesons and dou- 
bly heavy baryons. The ellipsis represents the kinetic energy terms which are 0{l/mQ) but 
do not explicitly break the heavy quark spin symmetry as well as higher order terms. The 
coefficient of the diquark chromomagnetic coupling differs from Ref. |l| by a factor of 1/2 
which will be explained below. The ground state heavy quark doublet consists of a spin-0 
meson, P, and spin-1 meson, P*, and the ground state doubly heavy baryon doublet consists 
of a spin-1/2 baryon, S, and spin-3/2 baryon, E*. Matrix elements of the spin-symmetry 
breaking operators in Eq. can be computed using the methods of Ref. [l| to make the 

following prediction for the hyperfine splittings: 

3 

m-£* — irtj] = -{nip* — nip) . (2) 

This .suit is also easily denved ,n the ,u.r, „,ode.. see e.g. Q. It can also be gene.ali.ed 

to baryons with more than one heavy flavor, but since there is no experimental evidence 
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for bcq baryons to date, we will restrict our discussion to the case of only one heavy quark 
flavor. 

Since the publication of Ref. 0] it has become clear that HQET is not the correct effective 
theory for hadrons containing two or more heavy quarks (or anti-quarks). This is because 
HQET is formulated as an expansion in Aqcd/'^Qj while a bound state containing two 
or more heavy quarks has two additional scales, rnqv, the typical momentum of the heavy 
quarks within the bound state, and rnqv"^, the typical kinetic energy of the heavy quarks. The 
appropriate theory for dealing with problems involving two heavy quarks is Non-Relativistic 
Quantum Chromodynamics (NRQCD) . This theory has the same Lagrangian as HQET 
but is formulated as an expansion in v rather than Aqcd/'^q- The goal of this paper is to 
re-derive the Lagrangian in Eq. (Q) and thereby the prediction in Eq. (0) using NRQCD 
as a starting point. We will also discuss higher order corrections to doubly heavy baryon 
hyperflne splittings. 

When Ref. Q was published the result in Eq. (j2)) was of academic interest because no 
doubly heavy baryons had been observed. Recently the SELEX experiment has published 



evidence for a number of ecu and ccd states p, M,l2|- One may worry about the applicability 
of our results to doubly charm baryons since NRQCD assumes a hierarchy of scales uiq 3> 
rnqv ^ rnQv"^, while in bound states with two c quarks or a c and c quark these scales are not 
so widely separated: rric ~ 1.4 GeV, rricV ~ 750 MeV and rricv'^ ~ Aqcd ~ 400 MeV. There 
are cases where the f ^ expansion converges very poorly for charmonium, see for example 
the calculations of Refs. P, 0], and we suspect this will also happen for some observables 
involving doubly charmed baryons. In fact, perturbative predictions for production rates 
jioj l and weak decays 11 1 seriously disagree with SELEX data 0, 0,!?!. However, using the 
known D* and D masses in Eq. (j21) yields ms* — ^ 105 MeV which deviates from the 
experimental value of 78 MeV by only ~ 30%, which is the size of the expected corrections. 
Lattice calculations of doubly heavy baryon hyperflne splittings are also approximately 80 
MeV This makes it seem plausible that the v expansion could be effective for this 
observable. 

The main point of the paper is to understand how approximate heavy quark-diquark 
symmetry arises within the framework of NRQCD and to lay the ground work for computing 
systematic corrections to the leading order prediction in Eq. ^ which are clearly necessary 
for understanding the spectrum of doubly charm baryons. An important aspect of this paper 
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is our novel method for deriving the Potential NRQCD (pNRQCD) Lagrangian directly from 
the NRQCD Lagrangian of Luke, Manohar and Rothstein j^, which is sometimes referred to 
as vNRQCD. In the QQ sector, auxiliary color-singlet and color-octet fields are introduced 
via a Hubbard- Stratonovich transformation to eliminate four-quark interaction terms from 
the vNRQCD Lagrangian. The quark and antiquark of vNRQCD are integrated out of the 
theory, leaving an effective action for the color-singlet and color-octet fields and usoft gluons 



which corresponds to the leading action of pNRQCD [1^. This helps clarify the relationship 
between the two formulations of NRQCD at the level of the classical Lagrangian. 

In section I, we rederive the pNRQCD action in the QQ sector to 0{v) using our for- 
malism. In section II, we extend this formalism to the QQ sector and derive the leading 
order spin symmetry breaking couplings of color antitriplet and sextet diquarks to usoft 
chromomagnetic fields. Section III contains a discussion of soft gluons and how are these 
are handled in our formalism. We also discuss higher order corrections to Eq. Conclu- 
sions are given in Section IV. While this work was in preparation, a pNRQCD treatment 
of baryons with two and three heavy quarks appeared in Ref. Their prediction for 

the hyperfine splitting agrees with Eq. 0. They also consider diquarks with two flavors of 
heavy quark and give the pNRQCD Lagrangian to O(f^). 

I. PNRQCD FROM VNRQCD 

In this section, we describe our method for obtaining an effective action for composite 
fields from the vNRQCD Lagrangian by introducing auxiliary fields and integrating out 
heavy quarks and antiquarks. We illustrate our method by rederiving the pNRQCD La- 



grangian 12], which is formulated in terms of composite color-singlet and color-octet fields. 

Our starting point is the formulation of NRQCD due to Luke, Manohar, Rothstein jj]. 
The effective theory contains several different fields: heavy quarks, whose kinetic energy and 
momentum are O(mQf^) and 0{mQv), respectively, soft quarks and gluons whose energy 
and momentum is 0{mQv) and ultra-soft (or usoft) quarks and gluons whose energy and 
momentum is 0{mQv'^). The scales uiq and mqv are removed from the effective theory using 
the trick of rephasing the NRQCD fields so that derivatives on the effective field theory fields 
bring only 0(mQf ^) momentum . Thus, soft fields and heavy quarks will carry a three 

vector label which indicates the 0{mQv) momentum carried by these fields, while usoft fields 
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carry no such label. We refer the reader to Ref. ^ for a detailed discussion of notation, 
power counting, and the renormalization group in this theory. 

Any light degrees of freedom in a quarkonium or a doubly heavy baryon have energy 
and momentum O(Aqcd)- For the physically irrelevant limit of extremely heavy quarks 
Aqcd ^ fngv^, while for charmonium or doubly charmed baryons we have Aqcd ~ fngv"^. 
In either case, the light degrees of freedom are described by usoft fields. Soft degrees of 
freedom appear only in loops when computing matrix elements between heavy quarkonia 
or doubly heavy baryon states. Therefore, it is possible to integrate out the soft loops and 
derive an effective action for only usoft and heavy degrees of freedom. This is essentially 
the philosophy of (weakly coupled) pNRQCD 

The validity of effective field theory containing only usoft and heavy degrees of freedom 
is not so clear in the context of the velocity renormalization group introduced in Ref. In 
this method for summing logs of v, the renormalization scale of soft loops and usoft loops 
are correlated: /is = rnqv"^ and /ius = ^qi^-, and v is evolved from v = 1 to v = v. As argued 
in Ref. both soft loops and usoft loops must be included throughout the evolution to 
correctly resum logarithms of v. Thus in vNRQCD, one must integrate out soft gluons at 
the scale mgv^, after summing logarithms of v. However, in the pNRQCD approach to 



summing logarithms 



12j |. soft gluons are integrated out at the intermediate scale rnqv. We 



will not be concerned with resumming logs of v and so this issue should not affect the results 
of this paper. 

The leading order vNRQCD Lagrangian in the QQ sector is 

- E T^4^<T"t-pxUT'x-p + ... , (3) 

where the ellipsis represents higher order corrections and terms including soft gluons. We use 
a Fierz transformation to project the potential onto color-singlet and color-octet channels 
and obtain position space potentials by Fourier tranforming with respect to the labels: 

= f dh^e'^P-''^' . (4) 
(p-q)2 J Airr 
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The Lagrangian is then 



where 



V^^{r) = -Cf^, V^'\r) =[^-CFj^. (6) 

Next we introduce auxihary fields into the Lagrangian via a Hubbard- Stratonovich trans- 
formation. We add to the Lagrangian 



A£ = 1 dhV^'\r) |^5t-^e-''--^^t^lqj _ ^ e'P-^-^ X-p V^p 

+ V^'\r) ||0r"t_^e-^q ry2^tTa£^j |^0r" - ^ e^P- X-p T'^ ^p ^ • (7) 

Note that the variable r is a continuous variable rather than a label. The terms quartic in 
the quark and antiquark fields cancel by construction and we are left with 



+ fd' rV^'\r) (slSr - Si e^^'^'-^X-p^p - S^J^ ^"^''^T^^S^^ ) 

V pVC qVcy 

+ jd'' rV^® (r) (^O^^O^ - Of ^ e*P"V2x-pT>p - ^ g-^'^ '-v^^jT^xlq j ■ 

Integrating out the fields Sr and yields the original NRQCD Lagrangian. Integrating out 
the fields ipp and Xp will yield an effective action for the fields and which can be used 
to study quarkonium at low energies. The terms in the Lagrangian of Eq. (jH)) which involve 
the quark and antiquark fields are reparametrization invariant [J, Q]. The terms involving 
the auxiliary fields are not, but can be made reparametrization invariant by making the 
following substitution: 



e^P-^X-p(T'^)^P ^ e'"-^ (x-pe-^-°/') {T'^) (e^'^/^^p) . (9) 

6 



In this expression both covariant derivatives act only on the fields inside their respective 
parantheses and in the first set of parentheses the covariant derivative acts on x-p from the 
right. In NRQCD power counting, p ~ rngv so we must take r ~ {mQv)~^. As a result 
r ■ D ~ f and the factors e^""'-^/^ in the Lagrangian can be expanded in powers of r ■ D. The 
0{v) term in this expansion includes a coupling to usoft gluon fields. When we integrate 
out the ipp and Xp fields, this coupling is needed to derive the couplings of the composite 
fields to usoft chromoelectric fields at 0{v). 

In passing we note that we can rewrite the exponentials of covariant derivatives using 
Wilson fines. We define 

Wr(x) = Pexp dsr-A(x + sr)^ , (10) 

which obviously obeys r • D Wr{x) — 0. Then we have 

r • DW^r(x)/(x) = VFr(x)r • a/(x) , (11) 

which is equivalent to the operator relation 

r • D = Wr(x) r • d Wj"(x) , (12) 

where r • D and r • d act only to the right. This allows us to write 

(e'-°/Vp(x)) = W-r(x) {e''-'/'W}{^)M^)) 

= iyr(x)H/;(x + r/2)^p(x + r/2) 

= Wr{x, X + r/2)V'p(x + r/2) , (13) 

where 

Wr{^,^ + T/2) = Pexp(^zgsJ^ ' ^ ("^ + ^ " 4)) ' ^^^^ 

is a Wilson line connecting the points x and x + r/2. Then the equations of motion for 
S'r(x) and 0"(x) yield 

5r(x) = 5^^e^P-^X-p(x-r/2)iyr(x-r/2,x + r/2)V'p(x + r/2) 

0»(x) = ^e'P""x-p(x-r/2)iyr(x-r/2,x) V2r"iyr(x,x + r/2)V'p(x + r/2). 
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FIG. 1: Feynman rules for the coupling of Sr and to quarks and antiquarks 

We see that once the corrections which restore reparametrization invariance are included 
to all orders the fields S'r(x) are nonlocal products of quark-antiquark fields separated by 
a distance r with center of mass coordinate x. The Wilson lines ensure that the nonlocal 
fields transform correctly under usoft gauge transformations, namely as 

Ipp^Vi:^ Sr^Sr Or^^V^'O^ (15) 

where V is in the fundamental representation and V is in the adjoint with V^T^'V = V^^tK 
The composite fields are very similar to interpolating fields introduced in Ref. for the 
construction of pNRQCD. Because the separation r is a short distance scale compared to 
the usoft scale we can perform a multipole expansion in r. 

We will first derive the kinetic terms for the and fields, for which we can neglect 
usoft gauge fields. The Feynman rules for the couplings of these composite fields to quarks 
and antiquarks are shown in Fig. ^ In these diagrams p and — p are the label momentum 
of the quark and antiquark, respectively. The vertex is nonvanishing only when the labels 
of the quark and antiquark fields are equal and opposite. An effective action for the Sj. 
and fields is obtained by integrating out the fields ipp and Xp? which amounts to doing 
one-loop graphs with quark anti-quark fields on the internal lines and an arbitrary number 
of external auxiliary field lines. This is similar to the theory of the superconductivity in 
which one introduces a pairin g fie ld and then integrates out the fermions to get an effective 
potential for the gap, see Ref. |l7|. However, in this case we are only interested in processes 
involving a single QQ pair and therefore only need the part of the action that is quadratic 
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in the composite fields. Therefore, computation of the effective action requires a one loop 
diagram with two external auxiliary fields, S'^S or 0"'^0"'. (We will find terms with 5*^0'* or 
O^^S* once we include external usoft gluons.) 




FIG. 2: One loop contributions to effective action for Sr and fields. 
The one loop graph in Fig. |21 with two external composite fields, Sr and Sr', yields 



^31 pil-(r-r') 

-y^'\r)-^-—^V('\r') . (16) 



(27r)3 ' _rnQE-ie 

In evaluating this graph the sum over virtual labels and integral over virtual residual mo- 
mentum have been combined into a single integral. The one loop diagram with external 
composite fields and O^, shown in Fig.|21yields a similar result. The effective Lagrangian 
obtained from integrating out the quark and antiquark fields is given by 

C = / d'vd'v'Sl[V^'\r)5\v-r')+ / ^K«(r)— ^ v^'\r')] Sr> (17) 

J \ J (2vr)'^ V/niQ - E - it J 

r / r d'^\ pjl-(r-r') \ 

+ / rf^rrfSr'Of v^^\r)5\r - r') + / T^V^^'^r)— ^ V'^^\r') O^, . 

J \ J (27r)^ V/niQ - E -le J 

This Lagrangian is most easily interpreted using the following notation. We think of the 
fields ^r, as vectors on a Hilbert space spanned by eigenkets of r, so we write Sr = (r|5'). 
The fields also depend on x but this dependence will be suppressed for simplicity. We define 
operators V^(^'^) whose matrix elements are (r'|V^'^^'^'' |r) = V^^'^\r)5^{r — r'), momentum 
eigenkets, defined by (r|l) = e*' "" and a free Hamiltonian, Ho\i) = P/mQ|l). In this quantum 
mechanical notation, the action can be written as 

C = « - \>«-^^ 1>«|^) + (C^I^X^) - K(^)|0") , (18) 

E — Hq E — Hq 

where E = id/dt. Inserting complete sets of states 1 = / d^r\v){r\ = J in the 

appropriate places in Eq. (fTH|) allows one to reproduce Eq. |T7|). Using algebra, it is easy to 



show 

m - V^^V\4>) = {^\E-H-iE- H)-J—-iE - HM) (19) 

= m-^^iE-HM. 

E - Hq 

where (f) denotes either composite field, V is the relevant potential, and H = Hq + V . The 
equation of motion is clearly [E — H)\(j)) = or 

z9o + — - Vir)] (t)ir) = . (20) 
J 

The second term on the right-hand side of the first line in Eq. ()19|) clearly vanishes by virtue 
of the equation of motion and can be removed by a field redefinition. Then the action is 
given by 

C = {S\E-Hq-V^^^\S) + {0^\E-Hq-V^^^\0^) 

= J dhSl (^ido + ^ - V^^^r)^ Sr + j dhO^^ (^ido + ^ - V^^^r)^ . (21) 



121. 



This is just the leading order Lagrangian of pNRQCD 

Next we consider corrections to the effective action that involve one loop diagrams with 
external usoft gluons. One complication is that the couplings of the Aq gluon to the elemen- 
tary quark and antiquark fields that come from the covariant derivatives in the kinetic terms 
of Eq. Q are not suppressed by any powers of v. (Couplings to Ai gluons are v suppressed.) 
Naively, this suggests that one loop graphs with an arbitrary number of external Aq gluons 
must be summed to lowest order in f . Of course, one can go to = gauge and all these 
diagrams disappear. Therefore, the only corrections to the effective action that involve the 
Aq gluons to lowest order in v are simply the terms required to make the derivative acting 
on Op in Eq. (PT|) covariant. In an arbitrary gauge there are cancellations between diagrams 
involving arbitrary gluons so that in the end the correction to the effective action is simple. 
The calculation of the effective action is greatly simplified if a field redefinition is performed 
that removes the 0(f°) interaction of the quark fields and the Aq gluon. This of course 
introduces new couplings to Aq gluons in higher order v suppressed interactions, but these 
can be treated perturbatively. 

Consider the kinetic terms for the quarks: 

^-E^i^'^o-^^^)"'.. (22) 
10 



To remove the Aq coupling we we make the field redefinition Q| 

= 5*0 = Pexp (^-igs j dsAQ{x + s)^ , (23) 

where 5*0 is a Wilson line obeying iDqSq = 0. Therefore we have SIiDqSq = ido and the 
kinetic term for the quarks in terms of the new fields is: 

^_^^,|^^^^_p^-apS|^-^j^-^ (24) 

We can write the terms with the Wilson lines more neatly if we use the following identities. 

^Jd.^o = + [^Jd.^o] = v. - zgB, = . (25) 
Brackets indicate that the covariant derivative acts only on terms inside the brackets. The 



the new fields the action is 



field Bj = ^[S'gDjS'o] is a gauge invariant generalization of the field Aj. Written in terms of 



^-i:^i(«o-^j4. (26) 

i.e., in terms of the new fields the Lagrangian is the same as before with the replacements 
—>■ , Aq —>■ , Ai Bj. When written in terms of the new fields the the action has no 
0(f°) interactions with gluons. The terms involving auxiliary fields are 



(27) 



+ j rV(«) (r) ^Of - Of ^ e*P" V2x-p3^o'^'^p ' OrJ2 e'*''" V2V^;3^o"^'X-q 

The terms involving S are the same as before once we make the substitutions if) ^ if) and x ~^ 
X- The terms with O" contain an additional octet Wilson line, defined by SqT°'Sq = y^^T^. 
This octet Wilson line can be absorbed by redefinition of the composite field: O^ = y^^'O"'. It 
is straightforward to show that after the field redefinition covariant derivatives D in Eq. 
are replaced with D. So the field redefinition can be implemented by making the following 
substitutions in the Lagrangian: 

^ V^, X ^ X, ^0 ^ 0, A, ^ B„ ^ S, O^'^O^ (28) 
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When rewritten in terms of the new fields, the action contains no f ° interactions involving 
gluons. We can now compute the effective action to 0{v) for the fields S, 0°^ and Bj. Once 
this is completed the action can be rewritten in terms of the original fields by undoing 
the field redefinitions if desired. The result is equivalent to what is obtained if one uses 
Aq = gauge and then imposes gauge invar iance. Performing the field redefinitions makes 
it possible to maintain gauge invariance at all stages of the calculation. 



k,a P,i 




FIG. 3: 0{v) couplings of composite fields to external gluons, quark and anti-quark fields. 

The relevant Feynman rules for graphs with one external usoft gluon, a quark, an anti- 
quark, and a composite Sy. or field are shown in Fig. |21 These are obtained by extending 
the Lagrangian in Eq. (jHJ to include the correction in Eq. Q which restores reparametriza- 
tion invariance and expanding to lowest order in the gluon fields. There are also 0{v) 
interactions in the coupling of quarks and antiquarks to usoft gluon fields. There are no 
other 0{v) interactions in the Lagrangian at this order. At higher orders in v, it is no longer 
sufficient to consider only the reparametrization invariant extension of the Coulomb inter- 
action in Eq. 0. At order O(f^) there are new potentials obtained from matching QCD 
onto pNRQCD Q- 

The four diagrams contributing to the 0{v) effective action for composite fields are shown 
in Fig. EJ Because of color conservation at least one of the external fields must be an octet. 
We will describe the evaluation of the graph with external Sr and O'^, in some detail. Let 
E and E' be the energy of the Sr and O'^, fields, respectively. Let the polarization index of 
the gluon be i and its color index b. Working in the rest frame of Sr we find that graph b). 
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a) 





c) J=^ 





FIG. 4: O(t') loop diagrams contributing to the coupling of composite fields and usoft gluon fields. 

in which the gluon connects directly to a quark line via the p ■ A coupling in Eq. (jH)), yields 
b) = f ^'"'^'^ ^'"'^''^ Fc-'-V-' (29) 

-Z V 16 6 



^ ^ =V^53(r-r'). 



The first line is obtained using the previously derived Feynman rules and performing the 
integral over the energy in the loop by contour integration. In the second line we have 
replaced factors of the momentum appearing in the integrand with derivatives of r and r' 
acting on the exponential factors in the integrand. In the last line we used the equations 
of motion to set E + '^l/niQ = V^^\r) and E' + 'Vl,/mQ = V^^\r'), then performed the 
integral over the loop momentum. Graph c), in which the gluon is attached to the antiquark 
line, gives the same result as graph b). Graph a) yields 

_ g r dH vW(^r)V('){r')e'^'-'^e^-^' ^^^^ 



^/2K J (2vr)3 E'-iymQ 

9 ^^.y f dH V-W(r)V-(^)(rO ^_,,,^,,,, 



^2Nl J (27r)3 E' + Vl/niQ 



9 rab 



where we have used the Feynman rules of Fig. IHl and the same manipulations used in the 
evaluation of graph b). Graph d) yields 
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V2N, 



■.6''^r"V^^\r')6''{r-r'). (31) 
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Combining all four diagrams then yields the following effective Lagrangian for the composite 
fields: 



r 



9 



—I 



V2K 

9 

V2K 



2 



2Vr' 



_r'\/(8)(r')) 



[VrOjl^r + r (y(i)(r) - \/(8)(r))Ot"5r 



(5=^(r - r') • B"Sr (32) 



In the second line we have integrated by parts and performed the integral over the 5-function. 
Now using the leading order equations of motion 

V(i)(r)Sr = (tdo + ^] Sr (33) 

and repeatedly integrating by parts one obtains 

= J= [ dhr- doB" 61" S, 



V2K 



(34) 



(35) 



To write this in terms of the original fields, we first note 

doB = -do[SlBSo] = -Sl[Do,-D]So = S^ESo, 
9 9 

where B and E are Lie algebra valued. In terms of components this equation can be written 
as (9oB" = — 3^o" E'', where 3^o" is the Wilson line defined shortly after Eq. pTjl . We can use 
this to rewrite Eq. (j35p in terms of the original fields using 3^Q''0^'*(r) = O^'^(r). 

Graphs of Fig. |3] with two external octet fields, O"^ and and an external gluon with 
color index c are the same as graphs with one external singlet field except for the replacement 
of the color factor 5^^/^/2Nl d"*72. The 0{y) effective Lagrangian for the composite 
fields is 

^- = 9^J dhTT[Olr-E]S, + h.c.+g J dhTrpUr ■E,Or}] . (36) 
Here we have defined = O'^T"^ and E = E"T" . This result is in agreement with Ref . 12 1 . 



II. AN EFFECTIVE ACTION FOR COMPOSITE DIQUARK FIELDS 



The NRQCD Lagrangian relevant to constructing an effective action for diquarks is 



(p-iD)- g 

+ o- ■ B Xp 



2mr 



2mc 



9l 



2 " IP - qj 

p,q 



Xq-^ XpX-q-^ X-p ~^ 



(37) 
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The ellipsis represents higher order corrections as well as terms including soft gluons. Also 
included is the leading spin symmetry violating interaction of the heavy antiquark since it is 
necesssary for generating hyperfine splittings. The following spin and color Fierz identities 

TtiTjk = ~o ^(^mij(^mik + 3 d^^'^^dl'P'^'^ , (38) 

m (mn) 

are used to the project the potential onto the 3 and 6 channels and to decompose the quark 
bilinears such as Xp X-p iiito operators with definite spin. In the first line, the Greek letters 
denote spin indices, the a* denotes the Pauli matrices and e = icr"^. In the second line, 
Roman letters refer to the color indices and the matrices d-™""^ are symmetric matrices in 
color space defined by 

^(rn.) ^ I + ^r^D/v^ rn^n ^ ^^^^ 

I ^r^j m = n 

In denoting the matrices by {mn) we do not distinguish {mn) from {nm) so there are six 
distinct matrices. After Fourier tranforming with respect to the labels the Lagrangian can 
be rewritten as 

C ^ >^^.^ + Ext(^^-^^ + ^'^-B)xp (40) 

- i| d'rV^'\r) (^e'^^-e.,,l{xi),ae{xl^),^ ■ (Y.e^''''leaM-^)iMXp)r 




where 



y(3)(,)= V^'\r) = \^. (41) 

6 r 6 r 



In Eq. fl40j) . we have left color indices explicit but suppressed spin indices. Note that the 
Fierz in Eq. ()38|1 introduces four terms but two vanish due to Fermi statistics. Diquarks 
in the 3(6) representation must be in a spin-1 (spin-0) state. Now we want to introduce 
composite fields which annihilate and create heavy diquarks. We define these fields as 

p 

= 5^e^--ri=dlr)(x-p).6^(Xp),. 



p 
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The interaction terms can be made reparametrization invariant using a replacement anal- 
ogous to that in Eq. 0. Note that T* is a vector and Sr'"'^'' is a scalar, as required by 
Fermi statistics. In the QQ, sector, we introduced fields Sr and that are bi-spinors which 
could be further reduced into their spin-0 and spin-1 components. If the diquarks are com- 
posed of two different heavy antiquarks, Fermi statistics no longer correlates the color and 
spin quantum numbers. We would have to either use fields in reducible representations of 
heavy quark spin SU{2), as is usually done in pNRQCD [3, ll^i or else introduce addi- 
tional fields for spin-0 3 and spin-1 6 diquarks. The diquark fields enter the theory via the 
Hubbard-Stratonovich trick: 

\ P 



X 



- Y: e--^ rfr^(X-p).e^(Xp), j . (42) 

The Feynman rules coupling diquarks to two heavy antiquarks from the interaction terms in 
Eq. are given in Fig. El The Greek letters a and (3 in these diagrams are spin indices of 
the antiquarks. The Feynman rules for the 0{v) couplings of diquarks, two antiquarks and 
a gluon which come from the reparametrization invariant extension of Eq. ()42|) are given in 
Fig. El 

Diagrams analagous to those in Fig.Hlgive rise to electric dipole couplings of the diquarks 
but in our case these couplings vanish. Electric dipole couplings do not fiip the spin of the 
heavy antiquark and therefore cannot couple diquarks of different spin such as T* and 
^(^mn) diquarks are composed of two different flavors of antiquarks, then electric dipole 
transitions between states with different color but identical spin are possible 13|. Electric 
dipole moments of the Tl and Tii"^"'^ vanish when the two antiquarks in the diquark have 
equal mass |13|]. Ref. [l3| has computed the electric dipole couplings in the case of non- 
identical antiquarks and we have checked that our methods reproduce their electric dipole 
couplings in this case. 

The chromomagnetic coupling of the diquarks is obtained by graphs like b) and c) of 
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P,3,a 
-p,k,p 

V2 



limn) 




FIG. 5: Feynman rules for the coupling of the diquarks T* and Sr™'"^ to antiquaries. 



m,b P,J,a 

X 

-p,j,/3 



FIG. 6: 0{v) couplings of diquarks to gluons and quarks. 

Fig. m with the p ■ A coupling of the antiquarks replaced with the cr ■ B coupling in Eq. (|4m) . 
These diagrams are O(f^) and give rise to the lowest order spin symmetry violating cou- 
plings of the diquarks. Other 0{v^ diquark couplings that do not violate heavy quark spin 
symmetry can be found in Ref . 13[ . For the lowest order chromomagnetic coupling we find 

More chromomagnetic couplings are possible if one considers diquarks composed of two quark 
flavors. For instance, the Sr™""^ field does not have a chromomagentic moment because it 
has spin 0, but if there are two flavors of heavy quark then the spin-1 field in the 6 of color 
has a nonvanishing magnetic moment jl3| . 

This Lagrangian is not yet in the form of Eq. (Q). To get an effective action for local 
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diquark fields of the form of Eq. (P) we can formally expand the composite fields in the 
eigenfunctions of the leading order Hamiltonian. For example, let 

T; = 5^T^0„(r), (44) 

where the 0ri(r) are an orthonormal basis of energy eigenfunctions satisfying 

f-Yl + V^'\r)] 0„(r) = -5„0„(r) . (45) 

Then the Lagrangian for nonlocal composite field T* can be rewritten as a local Lagrangian 
with an infinite number of fields: 

Ct = f dhTl f^Do + — - V^^\r)] Tr - f dhtTl B x 

J \ ruQ J 2mQ J 

= J2 TU^^o + 5n)T„ - ^ Tt . B X T„ . (46) 
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(Here we have suppressed SU{3) indices.) This assumes that the spectrum of the two-body 
Hamiltonian in Eq. is discrete. This is not the case within perturbation theory but 
must be true nonperturbatively in a confining theory like QCD. There are many heavy 
diquark fields and each has a different residual momentum given by the binding energy of 
the diquark in the nonrelativistic potential in Eq. ()45|1 . If one is interested in the properties 
of the ground state doubly heavy baryon only the lowest lying diquark field is relevant. 
Diquark fields corresponding to excited states of the diquark have an excitation energy of 
O(mQf^) as do the low lying states of the Sr'""'' fields, whose constituents feel a repulsive 
potential at short distances and therefore must be heavier by O(mQf^) than the ground 
state of the diquark in the 3 representation. For rriQV^ ^ Aqcd one can integrate out these 
degrees of freedom and one is left with a diquark Lagrangian of the form of Eq. with an 
additional residual mass term due to the diquark binding energy. This residual mass term 
does not affect the leading hyperfine sphtting prediction of Eq. ^ at 0(l/mQ). 

III. SOFT GLUONS 

In this section we discuss the role of soft gluons in the effective theory for composite fields. 
In vNRQCD soft gluons appear at leading order in v, where the coupling to the quarks is 
given by 

^soft = -Y E i^l'[^',,A''^']U,uip,p',q,q')^p + {i^^x)- (47) 
p,p',q,'i' 
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Here UnvijX p', q-, q') is a complicated function of the momenta which can be found in 
Refs. |l9|. At leading order in v the soft coupling respects heavy quark spin symme- 
try, with symmetry breaking contributions to t/(p, p', g, q') suppressed by one power of v. 
There are also soft couplings with an anti-commutator structure, {y4^,y4^,}, whose leading 
contributions are suppressed by one power of v and preserve heavy quark spin symmetry. 

Integrating out the heavy quarks via the Hubbard- Stratonovich transformation results 
in a coupling of the composite fields to soft gluons. However, soft gluons carry energy 
and momentum of 0{mQv) so emission of soft gluons couples diquark fields with different 
momentum labels. Because the effective theory is applied to diquarks which are bound 
in a baryon strong interactions can only change the composite field's total momentum by 
O(mQf^) or smaller. Thus soft gluons cannot appear as external degrees of freedom. They 
do, however, appear in loops where they can form a tadpole, give a higher order (in asimqv)) 
correction to the potential, or connect to an external usoft spectator. This last contribution 
is important for doubly heavy baryons since there is a valence usoft quark, and couplings 
of the diquark to the usoft quark via soft gluons could in principle lead to spin symmetry 
breaking interactions. 

One can try to derive a coupling of soft gluons to composite fields by evaluating a diagram 
like that in Fig. [7J where the gluons are soft rather than usoft. However, such a diagram is 
problematic because the soft gluon coupling depends on the label momenta, and this coupling 
becomes singular in the momentum region of interest. To be specific consider Fig. [7| Let 
one gluon carry incoming label three-momentum q and the other outgoing three-momentum 
q', and let the heavy quark and anti-quark coupling to the initial composite field have three 
momentum 1 and —1, respectively. The heavy quark coupling to the final state composite 
field has label momentum 1 + q — q'. Since the final heavy quark and antiquark have total 
nonvanishing label momentum q — q' the coupling to the final state composite fields vanishes 
except at q — q' = 0. When q — q' = the soft coupling is singular since U^tj{p,p',q,q') 
has terms proportional to q*)/|p — p'P, and |q— q'| = |p — p'| by label momentum 
conservation. 

The singular coupling arises because one is integrating out a gluon exchanged in the t- 
channel between the heavy quark and soft gluon fields appearing in Eq. ()47|). For generic 
p and p' this gluon has energy O(mQf^) and momentum 0{mQv). This mode is off-shell 
and should be integrated out of the theory. But when p = p', this mode has energy and 
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FIG. 7: Coupling a composite field to soft gluons. 

momentum O(mgf^) and is therefore a usoft mode which should not be integrated out of 
the theory. Therefore, the singular interaction can be removed from the theory by adopting 
the prescription that p = p' (and q = q') is excluded from the sum in Eq. (jTTj) . Then there 
will be no coupling of soft gluons to diquarks. 

Suppose as an alternative we adopt an ad-hoc prescription for dealing with the singularity, 
such as replacing (p— p')^ with (p— p')^-|-A^, where is an IR regulator, in the denominator 
of singular terms in Ufj,i^{p, p', q, q'). Then there will be a coupling of diquarks to soft gluons. 
However, the loop graphs with this coupling will vanish because the integrals over soft 
momentum are scaleless, since there is no external soft momentum, and scaleless integrals 
are set to zero in dimensional regularization. 

Integrating out the large seale ,n, with externa, heavy and nsoft quarks give rises to four 
quark operators of the form |20| 

+ ^^V^V + ^^V^^crV^ ■ 9 775? • (48) 

The operators with coefficients C2 and C4 violate heavy quark spin symmetry and can induce 
heavy quark spin symmetry violating couplings of the diquarks in the effective theory. Using 
the methods developed in this paper it is is straightforward to derive these couplings. For 
example we obtain 

C = ^ [dhtTl^-q j^.T'^q ^ x . (49) 

Such operators are suppressed by 0{as{mQ)'^v^) relative to the leading hyperfine splitting 
term. The matching coefficients, q, start at O(a^) and are calculated at the scale mq. 
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The usoft quark bilinear, q-y^^q, is 0{v^) while usoft chromomagnetic field, B, is 0{v'^) in 
vNRQCD power counting. Finally, there are corrections to the diquark chromomagnetic 
coupling coming from vNRQCD operators coupling four heavy quark fields and a usoft 

chromomagnetic field. The coefficients of these operators are suppressed by two powers of 
as and should give the leading correction to the prediction for the hyperfine splittings of 
doubly heavy baryons. 

IV. CONCLUSIONS 

Baryons with two heavy antiquarks can be understood as a point-like diquark, composed 
of the two heavy antiquarks, interacting with the light antiquark. The reason is that the 
heavy quarks form a nonrelativistic bound state of size r ~ l/(mQf), which is much smaller 
than 1/Aqcd the wavelength of the light degrees of freedom. As a consequence there is 
a heavy diquark-quark symmetry which relates doubly heavy baryons to heavy mesons. 
However, the dynamics of the doubly heavy baryons is described by NRQCD not HQET. 
In this paper we show how heavy diquark-quark symmetry arises within the framework 
of NRQCD. Starting from the QQ sector of the vNRQCD Lagrangian we use a Hubbard- 
Stratanovich transformation to derive a Lagrangian in terms of color-singlet and color-octet 
degrees of freedom. The Lagrangian we derive corresponds to the leading pNRQCD action. 
This analysis clarifies the relationship between vNRQCD and pNRQCD. We extend the 
formalism to the QQ sector, where we recast vNRQCD into a theory of diquarks in the 3 
and 6 representation of color. We include in our Lagrangian the leading symmetry breaking 
terms, which involve the chromomagnetic fields. In our derivation we consider the effects 
of soft gluons in vNRQCD, and argue that they only give corrections to potentials. In 
particular we claim that soft gluons do not contribute to the coupling of spectator quarks to 
the diquarks since any possible contributions vanishes in dimensional regularization. As a 
consequence operators coupling light quarks to diquarks only come from integrating out the 
scale niQ, and are suppressed by asimQ^v^ relative to the leading spin symmetry breaking 
term. 

To arrive at an effective action for diquark fields that is of the form of the HQET La- 
grangian we expand the diquark field in eigenstates of the Hamiltonian. When rnqv'^ 3> 
Aqcd we are justified in keeping only the lowest energy mode, and thereby obtain an HQET 
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like Lagrangian with an additional residual mass term due to the diquark binding energy. 
This additional term does not affect the leading hyperfine splitting prediction. Thus we are 
able to derive spin symmetry relations between heavy mesons and double heavy baryons 
from NRQCD. These spin symmetry relations are satisified at the 30% level for the charm 
sector even though the required hierarchy of scales rriQ ^ rrigv ^ rnqv"^ ^ Aqcd is not 
well satisfied for the cc system. The deviations from the spin symmetry relations are of the 
expected size. This makes a v expansion for this observable reasonable, and motivates a 
study of the higher order corrections which will be presented in a future publication. 
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